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The current-voltage (I-V) characteristics ol armchair graphene nanoribbons under a local uniaxial 
tension are investigated by using first principles quantum transport calculations. It is shown that 
for a given value of bias-voltage, the resulting current depends strongly on the applied tension. The 
observed trends are explained by means of changes in the band gaps of the nanoribbons due to the 
applied uniaxial tension. In the course of plastic deformation, the irreversible structural changes 
and derivation of carbon monatomic chains from graphene pieces can be monitored by two-probe 
transport measurements. 
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I. INTRODUCTION 

Graphene, as a (2D) monolayer honeycomb structure 
of carbon, has attracted a great deal of interest since its 
successful preparation in 2004.'^' Due to its unique me- 
chanical, structural and electronic properties, graphene 
have been realised as an important material for nu- 
merous theoretical investigations and promising applica- 
tions. Among these are charge carriers behaving as mass- 
less Dirac fermions,^ Klein tunneling,^ ^ ballistic trans- 
port at room temperature,^ ^ and anomalous quantum 
Hall effects.'^ From experimental points of view, field- 
effect transistorJ^, micromechanical resonatorJ^^, gas 
sensorJiS of graphene have already been proposed. Most 
of these are directly related with its transport properties. 

Earlier transport studies predict that spin-valve de- 
vices based on graphene nanoribbons can exhibit mag- 
netoresistance values that are thousands of times higher 
than previously reported experimental values .'i^ Unusual 
effects of dopings on the transport properties of graphene 
nanoribbons were also reported.fiSHlIl Nevertheless, the 
transport properties of graphene nanoribbons under uni- 
axial tension have not been fully explored even from the 
theoretical points of view. While the effect of strain 
on the electronic properties of graphene is becoming an 
active field of studyjl^^ the transport properties and I- 
V characteristics of nanoribbons under local or uniform 
strain is of crucial interest for development of future de- 
vice applications. 

In this study, based on state-of-the-art first-principles 
quantum transport calculations, we investigate the ef- 
fects of uniaxial strain on the current- voltage (I-V) char- 
acteristics of graphene nanoribbons. We showed that 
elastic strain can alter the electron transport properties 
dramatically. In some cases, under a 10% strain, the 
current can change as much as 400-500 %. However, the 
variation of current with strain is sample specific. Even 
more remarkable is that the chain fo rmati on of carbon 
atoms from the graphene nanoribbon d^^ l ^^ undergoing a 



plastic deformation can be monitored through I-V char- 
acteristics showing negative differential resistance. 

II. MODEL AND METHODOLOGY 

Geometry relaxations and electronic structures are 
calculated by using SIESTA package,ll2l which uses nu- 
merical atomic orbitals as basis sets and Troullier- 
Martin typ^^S! norm- conserving pseudopotentials. The 
exchange-correlation functional of the generalized gradi- 
ent approximation is represented by the Perdew-Burke- 
Ernzerhof approximation."^ A 300 Ryd mesh cut-off is 
chosen and the self-consistent calculations are performed 
with a mixing rate of 0.1. The convergence criterion for 
the density matrix is taken as 10"'*. Brillouin zone (BZ) 
sampling of the calculations have been determined after 
extensive convergence analysis. The conjugate gradient 
method is used to relax all the atoms until the maximum 
absolute force was less than 0.05 (eV/A). Interactions 
between adjacent graphene layers is hindered by a large 
spacing of ~10 A. 

The electronic transport properties are studied by the 
non-equilibrium Green's function (NEGF) techniques, 
within the Keldysh formalism, based on density func- 
tional theory (DFT) as implemented in the TranSIESTA 
(Ref. 112) module within the SIESTA (Ref. [HI pack- 
age. A single C-plus-polarization basis set is used. Test 
calculations with larger basis set and mesh cut-off were 
also performed, which give almost identical results. The 
current through the contact region was calculated using 
Landauer-Buttiker formulajl^Sl 

/(H) = Go T{E, Vb)dE, (1) 

where Gq — 2(e^//ij is the unit of quantum conductance 
and T{E, Vb) is the transmission probability of electrons 
incident at an energy E through the device under the po- 
tential bias Vb- The electrochemical potential difference 
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between the left and right electrodes is eVb = fJ-L — fJ-R- 



III. I-V CHARACTERISTICS UNDER ELASTIC 
STRAIN 

The band gaps of armchair graphene nanoribbons 
(AGNRs), which we consider in this study, depend on 
their widthspl which are conventionally specified accord- 
ing to the number of dimer lines, N in their primitive 
unit cell. AGNR(A^)'s are grouped into three families, 
namely N = 3m — 1 family having smallest band gaps, 
the N = 3™ family having medium band gaps, and 
N = 3m + 1 family having largest gaps, where m is an in- 
teger. Band gaps of each family decrease with increasing 
m and eventually goes to zero as m — >■ oo. AGNRs are 
nonmagnetic direct band gap semiconductors, which can, 
however, be modified by vacancied^ and impuritiesll^Sl 

The NEGF technique used to study the electronic 
transport employs a two-probe system; semi-infinite left- 
and right-electrode regions are in contact with a confined 
central scattering region. A two-probe system, specific 
to AGNR with = 8, but representative of any N, is 
shown in Fig. [T] (a). Both electrodes and the central 
region are made from AGNR(8) . Periodic boundary con- 
ditions were imposed on the plane perpendicular to the 
axis of the nanoribbon. The carbon atoms at the edges 
are saturated with H atoms. The central region con- 
tains 5 primitive unit cells, with a total length of 21.76 
A (= 5co). The length of the central region is sufficient 
enough to avoid an abrupt change in electronic struc- 
ture while progressing from the electrode region to the 
strained region of interest. 

We first consider the electronic transport properties of 
the unstrained two-probe system presented in Fig. [T] (a) . 
To provide for an intuitive understanding of the trans- 
port phenomena, the band structure of the electrodes 
or the scattering region in their primitive unit cell are 
shown in Fig. [T] (b) . The lowest conduction and highest 
valance bands originate from tt*- and 7r-states, respec- 
tively. Unlike the perfect 2D graphene, where tt- and 
7r*-bands cross at the iC-corners of BZ, AGNR(8) is a 
direct band gap material having 0.20 eV band gap value. 
The calculated zero-bias transmission spectrum is given 
in Fig. nl (c), which apparently mimics the band struc- 
ture of AGNR(8). There is a region of zero-transmission 
with a width of 0.20 eV and located around the Fermi 
level, coinciding with the band gap of AGNR(8). Like- 
wise, the step-like behavior of the spectrum is related 
with the available conductance channels due to bands. 
The current as a function of the applied bias voltage Vb 
is presented in Fig. [l] (d). For this type of calculations, 
we increased VJ, in steps of 0.1 V and used the converged 
density matrix of the previous state as an initial guess for 
the next step. Applying a bias voltage shifts the Fermi 
level of the left-electrode with respect to the Fermi level 
of the right-electrode. The current starts flowing once 
the top of the valence band of the left-electrode matches 
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Figure 1: (Color Online) (a) Schematic view of two-probe 
armchair graphene nanoribbon system AGNR(8), having 
7V=8 dimers along the z-axis. Central scattering region, left 
and right electrodes are indicated. Carbons atoms are repre- 
sented by large (brown) and H atoms by small (light) balls. 
Primitive unit cell of electrodes and scattering region are the 
same and represented by dashed lines. The lattice constant of 
the primitive unit cell is co, and that of the central scattering 
region is 5co. (b) Band structure of AGNR(8) in its primitive 
unit cell, (c) Transmission spectrum of the system shown in 
(a) under zero-bias voltage, (d) I-V plot of AGNR(8) for a 
bias voltage from to 2.5 V. (e) The transmission spectrums 
for 4 different bias voltage calculated for the system shown in 
(a). 



in energy with the bottom of the conduction band of the 
right-electrode, as expected from the T(E, Vh — 0) for 
low-bias values. The T(E,Vb) does not alter much with 
the bias, since it is a uniform system and no significant 
permanent charge migrations should occur. This is evi- 
dent from the linear response of the current to the bias 
voltage for values of Vi, > 0.2 eV. As the calculation 
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of the current is very time consuming, the bias range is 
Umited from to 2.5 V. In Fig. [T] (e) the transmission 
spectrums for bias voltages of 0.6 V, 1.2 V, 1.8 V and 2.4 
V are also presented. As readily seen, the transmission 
T(E,Vf,) contributing to the current always keep near 1 
Go, and the higher values in transmission values move 
further away from the Fermi energy. This is due to the 
fact that as the bands of the leads move up or down in en- 
ergy scale with the varying bias, only a single conduction 
channel is open, or in other words, only one band crosses 
the energy of interest, at either one or both of the leads. 
This holds true for the bias voltages we consider in this 
study and as a result, we see a linear current response to 
voltage for zero-strain. 

Earlier, we have investigated the elastic and plastic de- 
formation of graphene and its nanoribbons under uniaxial 
tension. Mechanical properties were revealed from the 
strain energy, Es = ET{e) — Exit — 0); namely, the total 
energy at a given uniaxial strain e minus the total energy 
at zero-strain. Here, the uniaxial strain is e = Ac/cq, 
where cq and c = cq -\- Ac are equilibrium and stretched 
lattice constants of the nanoribbon, respectively. The 
tension force. Ft = —dEs{e)/dc and the force constant 
K — d^Eg/dc^ are obtained from the strain energy. Cal- 
culated elastic constants were in good agreement with 
available experimental data obtained from graphene.l^ 
Here we consider the I-V characteristics of AGNR(8) un- 
der a uniform uniaxial tension of the central scattering 
region for < e < 0.18. The strain is introduced as 
follows: The electrode atoms are fixed in their equilib- 
rium positions while the length of the central region is 
increased uniformly by Ac. Subsequently, the structure 
of the central region is fully optimized in a larger super- 
cell containing also unstrained electrode regions. In this 
respect, our study reveals the effect of a local strain in 
a long unstrained nanoribbon. The total energy of the 
system is recalculated. The strain energy Eg, is obtained 
according to above definition. Eg versus the elongation 
Ac, as well as e = Ac/5co plot for AGNR(8) system is 
given in Fig. [2] (a). The segment of AGNR(8) in the cen- 
tral scattering region undergoes an elastic deformation 
up to strain values e ~ 0.18, where the honeycomb-like 
structure is maintained, and the system returns to its 
original configuration if the tension is released. However, 
for higher values of strain, the system deforms plastically, 
where irreversible structural changes occur and the strain 
energy suddenly drops. Further information for this type 
of elastic and plastic deformation can be found in ReffTSl 

In Fig. [2] (b)-(d), we present the I-V plots of stretched 
nanoribbons. The electrode regions are identical to the 
no-strain case, but the central region under strain causes 
the changes. Once again, due to the band gap of the 
electrodes of 0.2 eV in AGNR(8), no current is observed 
up to a bias voltage of 0.2 V. The current response to 
bias voltage is linear for low strain, but becomes in- 
creasingly non-linear for higher strain. It is important 
to notice that higher strain in the central region induces 
stronger non-uniformity on its geometry and thus on its 
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Figure 2: (Color Online) (a) The strain energy Es, as a func- 
tion of elongation Ac, or strain e = Ac/co for AGNR(8) sys- 
tem, (b) I-V characteristics of AGNR(8) for a bias voltage 
from to 2.5 V for different values of strain. The increase 
and decrease trends are shown by arrows, (c) and (d) same 
as (b) for different systems, AGNR(6) and AGNR (10), re- 
spectively. 



electronic structure as compared to the electrodes. Equi- 
librium charge transfer may occur and alter the systems 
response to the non-equilibrium perturbation. This will 
result in a varying T(E,Vb) for different values of 14 . It 
is informative to compare the current values for systems 
under different strain at a given bias voltage. For exam- 
ple, at Vb=2 V, the current is around 62 fiA at e = 0. 
It increases to 65 fiA at e = 0.02, but steadily decreases 
for higher strain, having values of 28 at e = 0.12 and 
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Figure 3: (Color Online) The values of band gaps and currents 
of (a) AGNR(8) and (b) AGNR(IO) systems calculated as a 
function of strain, e. The currents are calculated for 2.0 V 
bias-voltage, (c) The band structures of AGNR(IO) under 
different strains. The band gaps are shaded and Fermi energy 
is set to zero. 



12 fiA at e = 0.18. One also notes that the I-V curve in 
Fig. [l] (d), which is almost linear for e = 0, starts to lose 
its linearity for higher values of strain as seen from Fig. [2] 

(b) . 

Other ribbons such as AGNR(6) and AGNR(IO) whose 
I-V characteristics are given in Fig. [2] (c) and (d). 
AGNR(6) belongs to the N = 3m family and it has a 
larger band gap (~1.04 eV) as compared to AGNR(8). 
As a result, we do not observe a current until 1 eV as 
seen from Fig. [2](c). The current steadily decreases un- 
til e — 0.10, then starts to increase as seen from Fig. [2] 

(c) . AGNR(IO) is another system which has a band gap 
value around 1.00 eV and its I-V characteristics are given 
in Fig. [2](d). In contrast to AGNR(6), the current first 
increases until e — 0.08 and then starts to decrease for 
higher strain values. All these results in Fig. [2] (b),(c), 
and (d) show that the current passing through nanorib- 
bons is very sensitive to the strain values and the behav- 
iors of I-V curves are sample specific. 

The increase and decrease of the currents given in 
Fig. [2] due to the changes in the strain is directly related 
with the electronic structure of the central scattering re- 
gion, which is modified as a result of changes in atomic 
structure under tension.'^^'^^' In Fig. [s] (a) we show the 
variation of current and band gap of AGNR(8) with ap- 
plied uniaxial strain. Here current values are extracted 
from Fig. [2] (b) for 2 V bias voltage. As seen from the 
plots, there is an inverse relationship between the cur- 
rent and band gap values. Any increase in the band gap 
decreases the current and vice versa. The same analy- 
sis performed for AGNR(IO) in Fig. [s] (b) also confirms 
this relationship. The changes in the band structures of 
AGNR(IO) can also be followed from Fig.|3](c), where the 



lowest conduction and highest valence bands approach to 
each other until e ~ 0.09 and then move away for higher 
values of strain. The band gap variations occur due to 
different nature of bands around the conduction and va- 
lence band edges exhibiting different shifts with strain. 
In particular, note that tt*- and 7r-bands of AGNR(IO) 
cross linearly at F-point by closing the band gap. This is 
the realization of massless Dirac Fermion behavior in a 
nanoribbon, which is semiconductor under zero-strain.^ 
In a simple model, an electron is ejected from the left 
electrode at an energy value lower than the shifted va- 
lence band maximum, for available ranges within the bias 
voltage. It is incident upon the central region, with lower 
chemical potential, and tunnels through to the right elec- 
trode, still lower in energy. The smaller the band gap 
value for the central region, the larger the number of 
possible states that participate in the tunneling, thus the 
larger is the value of the current. 



IV. TRANSPORT PROPERTIES OF AGNR(8) 
UNDER PLASTIC DEFORMATION 

While the elastic deformation imposes changes in the 
band gap and current values, the onset of plastic defor- 
mation results in dramatic changes in the structure. Af- 
ter yielding, the modification of honeycomb structure is 
somehow stochastic and sample specific. It depends on 
the conditions, such as the defects in the sample, the 
temperature effects and the rate of stretching. However, 
it has been shown theoretically— ^ that under certain cir- 
cumstances a long carbon atomic chair!^ (identified as 
cumulene having double bonds and polyyne with alter- 
nating triple and single bonds) can form in the course 
of plastic deformation of graphene, unless the edges of 
AGNR is not terminated with hydrogen. Upon further 
stretching, each carbon atom of graphene implemented 
to chain results in a stepwise elongation of the chain be- 
tween two graphene pieces. Monatomic carbon chain, 
which was derived experimentally from graphene,!^ can 
be a potential nanostructure for various future applica- 
tions. The important issue we address here is how these 
sequential structural changes reflects the transport prop- 
erties. 

In Fig. [4] (a) we present the atomic structure of a two- 
probe graphene nanoribbon system which is formed after 
the plastic deformation of AGNR(8) nanoribbon. A short 
chain containing 4 carbon atoms between the graphene 
flakes is formed in the scattering region and its zero-bias 
transmission spectrum is presented in Fig. |4] (b). This 
spectrum is composed of peaks rather than step-like lev- 
els as in Fig. [1] (c). The calculated I-V plot in Fig. [I] 
(c) also contains some peaks, which may lead to neg- 
ative differential resistance.ESl Similar situation also ex- 
ists for a longer chain in Fig. |4] (d), which occurred at 
a more advanced stage of plastic deformation whereby 
the nanoribbon in the central region is more stretched 
than in Fig. H (a). At the end, two more carbon atoms 
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Figure 4: (Color Online) (a) Schematic view of AGNR(8) system which is deformed plastically due to the high values of strain 
(e > 0.20). Its zero-bias transmission (b), and I-V characteristics (c), are given on the right. (d)-(e)-(f) are the same as 
(a)-(b)-(c) for a longer chain between graphene pieces. 



arc included to the chain. The differences between zero- 
bias transmission curves in Fig. and (e) occur be- 
cause of the energy level diagram and their positions rel- 
ative to Fermi levels are different. Also, the I-V curve 
corresponding to two carbon chains of different lengths 
occurring in subsequent stages of stretching are rather 
different. We note that the conductance of the longer 
carbon chain in Fig. and the corresponding current 
values (I) of a given bias voltage (V) can be higher than 
the shorter chain in Fig. Qa). This paradoxical situa- 
tion is related with the fact that some energies of the 
channels can be closer to the Fermi level as the length 
of the chain increases.!^ Further stretching of the sys- 
tem shown in Fig. can produce longer carbon chain 
structures. The length of these chains can be as long as 
10 carbon chains. As found for the structures in Fig.|4j^a) 
and Fig. |4|^d), the I-V characteristics of the longer car- 
bon chains will be different and will allow one to monitor 
the structural changes. Finally the plastic deformation 
terminates upon breaking of the chain. 



probe system can be modified with uniaxial strain. The 
current under a fixed bias can change several times with 
applied uniaxial strain. However, these changes are sam- 
ple specific and related with strain induced changes in the 
electronic structure near the band gap. Irreversible struc- 
tural changes and the formation of monatomic carbon 
chain between graphene pieces in the advanced stages of 
plastic deformation can be monitored through two-probe 
transport experiments. We believe that our findings are 
of crucial importance for recent active studies aiming to 
reveal the effects of strain on the electronic properties of 
graphene. Also our results suggest that these systems 
can be used as nanoscale strain gauge devices. 
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